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Integrable correlated electron systems (see, e.g. |J) with open boundary conditions constitute an important class 
of physically significant models which are exactly solvable by means of the Bethe ansatz technique and the boundary 
quantum inverse scattering method (QISM) [^,^). In this letter, we determine integrable boundary conditions for 
the g-deformed (or two-parameter) supersymmetric U model |6| of strongly correlated electrons. We present four 
classes of boundary conditions for this model, all of which are integrable by the supersymmetric boundary QISM 
formulated recently in Q . We solve the boundary systems using the coordinate Bethe ansatz method and derive the 
Bethe ansatz equations for all four cases. 

Let Cj ;(T and cj g denote fermionic creation and annihilation operators with spin a at site j, which satisfy the 

anti-commutation relations given by {c\ a , Cj, T } — 6ij6 aT , where i,j = 1,2, • • • ,L and a, r =|, |. We consider open- 
boundary g-deformed supersymmetric U models with the Hamiltonian of the form: 

L-l 
3=1 

where Bu (B rt ) stands for left (right) boundary term whose explicit forms are spelled out below, and Hjj +1 is the 
bulk Hamiltonian density of the g-deformed supersymmetric U model [^,^) 

H fj+i = - Y^ c \° c i+ la + ex P(^^(^ ~ si g n ( cr )7)%- -<r - \{n + sign(cr)7)n i +i,_ (7 ) 
u , 

+tp(c] T c]_ l Cj + i4.Cj + if + h.c.) + e~ 2a rij + e 2a n J+1 . (2) 

Here sign(cr) = 1 for a =f and —1 for a =j, iij a is the density operator rij a = c^ a Cj a , rij = njf + riji, and 77, 7 are 
two free parameters; and 



t p = — = [2e v (cosh j] — cosh 7)]" 



1 [ln[ Sinh ^ 7? + 7 ^ 1 - 7} (3) 
4 sinh i(ry — 7) 

Let us point out that in the limit of r\ — 7 (first discarding the chemical potential terms and then taking the limit) , 
the bulk model Hamiltonian (||) reduces to that of Bariev Q. This implies that, up to the chemical potential terms, 
the bulk model (||) is a generalization of the Bariev's model . 

We propose the following four classes of boundary conditions: 
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Case A : Bu 
B r t 

Case B : Bu 
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B r t 

Case D : Bu 
B rf 
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sinh 7(2 ^ +) \ sinh 4 
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sinh 2a / 
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sinh 2a 



sinh(7C 4 



(4) 
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sinh 7 ^ 2 2 ^ ^ ysinh „ 



sinh7 _ T n_i^i 
— ni T n u - e 7l > 2 ; ni 



sinh 2a 



sinh(7C + — a) 
sinh 2a 



e 7C +- a n Li ; 



e 7C - +a n U , 



sinh(7C_ + a) 

sinh 2a / sinh 7 



sinh 2^2^5+2 ^sinh^i 



(0) 



(7) 



where £± and (± are some parameters describing the boundary effects. 

Quantum integrabiiity of the system defined by Hamiltonian ([!]) with all four boundary conditions can be established 
as follows by means of the supersymmetric boundary QISM || in which the above Case A is treated in details as an 
example. We first search for boundary K-matrices which satisfy the graded reflection equations Q : 

1 22 1 

i?i 2 (ui - u 2 ) K- (ui)R 2 i(ui + u 2 ) K- (u 2 ) =K- (u 2 )Ri 2 (ui + u 2 ) K- (ui)R 2 i(ui - it 2 ), 



1 l 12 

Ri2{-U! + u 2 ) K + ( Ul ) M R 2 x(-u 1 -u 2 + 2) MK+ (u 2 ) 



12 1 1 

--MK+ (u 2 )R 12 (- Ul -u 2 + 2) M~ 1 K+ {u 1 )R 2 i(-u 1 + u 2 ) 



(8) 



where R(u) is the R-matrix of the g-deformed supersymmetric U model ||, R 2 \{u) — P\ 2 R\ 2 (u)Pi 2 with P being the 
graded permutation operator and M = diag(l, 1, e 27 , e 27 ) is the so-called crossing matrix. There are two different 
diagonal boundary K-matrices, KL(u), K^{u), which solve the first reflection equation in (||) |J: 



Ki(u) = 



sinh sinh 



7(g-~2) 



^diag {A- (u) , B_ (u) , B_ (u) , C_ (u) ) : 



K H (u) = diag(A_(u),A_(u),B_(u),B_(u)) , 



(9) 
(10) 



where 



B_{u 
C-(u 
A_(u 
B_(u 

A (II) 



sinh 



7(£- + u) ginh 7("~2 + g-) 



sinh sinh 2(!tz|±iO i 

e- 7 "sinh^-" M) sinh^- 2 + ^ ) . 



e -7« _ 1 



2 sinh(7C- + a) 
e 7 " - 1 



o7C-+a 



,7C-+a 



(11) 



2 sinh(7C- + a) 

The corresponding K-matrices K 1 ^ 11 ' ' (u) which obey the second reflection equation in (JsJ) are derived by isomorphism: 

K' + {u) = MKi(-u + 1), K^(u) = MK H (-u + 1) (12) 
We form the boundary transfer matrix t(u): 
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t(u) = str[K + (u)T_(u)K-(u)TZ 1 (-u)], T_(u) = R 0L (u) ■ ■ ■ R Q i(u). (13) 

Since K±(u) can be taken as K±(u) or K±{u), respectively, we have four possible choices of boundary transfer 
matrices: 

t A {u) = s tr[4( u )T_( u )i{[(i 1 )r 1 (-«)], 

t B (u) = str[K^ {u)T_{u)K n '(ttjTr^-u)], 
t c {u) = str[K% (u)T-(u)Kl{u)Tz\-u)\, 

t D ( u ) = str[Kl{u)T^{u)Kl\u)TI l {-u)l (14) 

which reflects the fact that the boundary conditions on the left end and on the right end of the open lattice chain 
are independent. Then it can be shown || that Hamiltonians corresponding to all four boundary conditions can be 
embedded into the above four boundary transfer matrices, respectively. We thus arrive at the four possible cases 
(H~@)i a ll °f which are compatible with the bulk integrability. 

We now solve the above boundary systems by using the coordinate space Bethe ansatz method. Following |u|f§, 
we assume that the eigenfunction of Hamiltonian ([[]) takes the form 

l*) = Yl 0=i. •• -,XN)cl iai •••4 jvctjv |o), 

N 

*<7i,-,c7 W (aJi, • • -,xn) = /^2^P A a Q1 ,---,a QN {kpQi, ■ ■ ■ ,kpQ N )exp(i^2k Pj Xj), (15) 

P 3=1 

where the summation is taken over all permutations and negations of ki , ■ ■ ■ , fcjv , and Q is the permutation of the N 
particles such that 1 < xqi < • • • < xqn < L. The symbol ep is a sign factor ±1 and changes its sign under each 
'mutation'. Substituting the wavefunction into the eigenvalue equation H^) = E\ty), one gets 

A. ■■ ,(7^ , • • • (' ' ' : kj , ki , • • •) — Sij {ki , kj )-A... )(7i5( 7, ... (• • • , ki, kj , • • •) , 
( — kj , ■ ■ ■) — s [kj] pi CFi )A a - ir ..(kj, ■ ■ •), 
A.. )CT( (- • • , -kj) = s R {k J ; PLat )A...^(- ■ ■ , kj), (16) 

where Sij(ki,kj) are the two-particle scattering matrices, 

S%j{ki : fcj)xi = Sijiki, kj)^ = I? 



S (k- k ) 12 - <? ■ (k- k) 21 - sm ( A » A j) 

sm(Ai — \j — i"f) 



sm 17 



sin(Ai - \j - £7) ' 

b v {k,,k 3 ) 12 -e sin(Aj — A 3 - — ij) Uj 

with Aj being suitable particle rapidities related to the quasi-momenta kj of the electrons by 

k(X) = 2 arctan(coth a tan A) (18) 
and s L (kj;pi ai ), s R (kj;pL ai ) are the boundary scatering matrices, 

1 - Piv, e' kj 



Afe j ;^)= 1 [ ^2^ e^(^) (19) 



with and pL CTi being given by the following formulae, corresponding to the four cases, respectively, 
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„ . 2a sinh2a _~fi_«-i 

Case A : p 1T = p u = pi = -e 2a e 7(1 ^, 

sinh - ^ - 

,„ sinh 2a , . 

p iT = P£1 = p L = - e - 2a - . -ETT e7( w) ; ( 2 °) 

sinh fv 3 ; 

Case B : p 1T = -e 2 ", p u = -e 2 " + . f^a e ^-+°, 

smh(7C- + a) 

Pit = -^ 2Q > Plx = -e- 2a - . , S ; nl ; 2a ^ (21) 

sinh (74+ — a) 

^ ^ 2a sinh 2a 

Case C : p 1T = p u = pi = -e 2a e 7(1 r), 

sinh - 2 ** ' 

m = - e - 2a , m - - e - 2a - . , s ; nl ; 2a ^ c +- a ; (22) 

smh(7C+ - a) 

Case D : p 1T = -e 2a , p u = -e 2 « + . S ^ 2 ° e*-+°, 

smh(7C- + a) 

,„ sinh 2a „n «+-, . . 

p LT =p Li =p £ ^-e- 2a - - ■ 7(2 _ g+) e 7(1 -T- ) . (23) 

2 



sinh 



As is seen above, the two-particle S-matrix (|l7|) is nothing but the R-matrix of the six-vertex model in the homoge- 
nous gradation and thus satisfies the quantum Yang-Baxter equation (QYBE), 

Sij(ki,kj)Sii(ki,ki)Sji(kj,ki) = Sji(kj , ki)Sn(ki, k{)Sij(ki, kj). (24) 

It can be checked that the boundary scattering matrices s L and s R obey the reflection equations: 

Sji{~kj, ki)s [kj)Piuj)8ij{, ki,kj)s (A^jpi^) 

= s {ki' 1 pi ai )Sji( kj,ki)s (kj' pi ai )Sij (ki, kj) , 
Sji ( kj, ki)s R (kj ; PL<jj )Sij {ki , kj)s (fcj ; PLoi ) 

= s R (k i ]p Lai )S : j l (k : j,-ki)s R (kj]p L(Jz )- 1 p ai )S : j l {k : j 1 ki). (25) 

This is seen as follows. One introduces the notation 

1 - pe~ tk 
1 — pe 1 



s(k;p)= - y _ lk ■ (26) 



Then the boundary scattering matrices s L (kj\p\ (Ti ) 1 s R {kj\pLe i ) can be written as, corresponding to the four cases, 
respectively, 

Case A : s L (kj]pi <Ji ) — s(—kj;pi)I, 

s R {k f ,p Lai ) = e^ L+1 ^s(k 3 ;p L )I; (27) 



c 2iXi sin(^-+Aj) q 

Case B : s L {kj]pi ai ) — s(— kj\p\{) | s^ll^^x^ 

* 1 

s R {k 3 -p Lai ) = e^ L+l h{k J -p Li ) [ 6 ' sin(2a+A 3 .) ) ; (28) 
Case C : s L (fc J ;pi (Ti ) = s(— kj',pi)I, 

s R (k j ;p Lai )=e i M L+ V S (k j ;p Ll ) ( 6 ' sin(2S± + A 3 .) ) ; (20) 

c 2iA,- sin( 7 i^+A J ) g 
Case D : S (feyJPlo*) = s(— kj-,p\{) \ sin( T -^-A i ) 

1 1 

S *(fc j; p iCTi ) = e lk ^ L+1 h{k f) p L )I. (30) 
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Here I stands for 2 x 2 identity matrix and p%i, pl[ are the ones given in ((2^). We immediately see that 
trivial solutions of the reflection equations (pq), whereas ( |28| ) are the diagonal solutions PJll[. 

The diagonalization of Hamiltonian (|l|) reduces to solving the following matrix eigenvalue equation 

Tjt = t, j = l,---,N, 

where t denotes an eigenvector on the space of the spin variables and Tj takes the form 

Ti = SJ (kj^i-kj ; pi*, )RJ (k,)R+ (kj^ikj ; PLaj )S+ (k 3 ) 

with 



are the 



(31) 



(32) 



s ?( k j) = S jtN (kj,k N ) ■ ■ ■ Sj,j+i(kj,kj+i), 
&j \ k j) = fej'-i) ' ' ' 

Rj i k j) = —kj) ■ ■ ■ Sj-ij(kj-i, —kj), 

Rt( k i) = Sj+ij(k j+ i,-kj) ■ ■ ■ S Nd (k N , -kj). 



(33) 



This problem can be solved using the algebraic Bethe ansatz method. The Bethe ansatz equations for all the four 

cases are 



Pla,PLa) 



|t sin^ ■ - A a - f ) sin(A 3 -+A a -f) 
W sin(A, -K + l i) sin(A J + A Q + ^) ' 



n 



sin(A a - Aj + f ) sin(A a + Xj + f ) 
!; sin(A« - Aj - ^) sin(A Q + A, - 



A/ 



= G(A a ;p la ,p LrT ) Yl 



13=1 

f3^a 



sin(A a - Afi + ij) sin(A a + A + ry) 
sin(A Q - A/3 - ij) sin(A Q + Ap - ij) ' 



(34) 



where 



F(kj\pi ai p La ) 



G(A a ;pi a ,p La ) 



s{kj\pi)s(kj-,p L ) (Case A) 

s(kj;pi^)s(kj;pLi) (Case B) 

s(kj]pi)s(kj;pL J [) (Case C) 

s(kj;pii)s(kj;p L ) (Case D) , 

1 

sin(A Q — ijC++i-%) sin(A a — 17^-+^) 
sin(A Q +i7^+— i^) sin(A a +?7^_ — 
sin(A et -t7C + +i 3f) 

sin(A„+i7C+— i4) 
sin(A Q — 27^- +2^-) 

sin(A Q 4-*7C- — ^j) 



(Case A) 
(Case B) 

(Case C) 

(Case D) , 



(35) 



where again pij , p^ arc the ones given in (|2lj). The energy eigenvalue E of the model is given by E = —2 52j=i cos k j 
(up to some additive constants, which we have dropped). 
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